Phase transition in a spring-block model of surface fracture 
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A simple and robust spring-block model obeying threshold dynamics is introduced to study surface 
fracture of an overlayer subject to stress induced by adhesion to a substrate. We find a novel phase 
transition in the crack morphology and fragment-size statistics when the strain and the substrate 
coupling are varied. Across the transition, the cracks display in succession short-range, power-law 
and long-range correlations. The study of stress release prior to cracking yields useful information 
on the cracking process. 
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Apart from its significance in industrial applications, 
cracking of a brittle material contains interesting physics 
and has attracted a lot of attention recently [Q. Two 
kinds of ongoing research can be identified: the first deals 
with the dynamic instability in the propagation of a sin- 
gle crack [|J, and the second is concerned with the col- 
lective behavior of many interacting cracks, such as the 
global pattern and distribution of fragments [||-|| . The 
latter receives less attention and will be our focus here. 

Although crack patterns in nature occur in a variety 
of contexts Q over a wide range of sizes from the mil- 
limeters on a monolayer of packed polystyrene spheres ||] 
to the kilometers of giant crack networks on playas, they 
are formed by the same basic processes. Typically, when 
an overlayer dries and shrinks, adhesion to a substrate 
resists shrinkage and induces internal stress. The stress 
may be released by slipping at the interface of contact 
or by cracking the overlayer; their competition leads to 
the variety of crack patterns. In recent experiments, the 
scale and geometry of patterns were found to depend on 
the strength of adhesion, the thickness of the layer and 
boundary conditions . Analytical approaches are dif- 
ficult and scarce, due to the complexity of multi-crack 
interactions j7|. Thus, simulations using simple models 
are expected to provide useful information p,H,B[. 

In this Letter, we present a spring-block model of frac- 
ture for a two-dimensional overlayer in contact with a 
substrate. Our goal is to identify the major control pa- 
rameters and to explore various statistical properties of 
fracture not readily accessible to experiments. A novel 
phase transition is found, whose origin can be traced to 
the interplay between stress redistributions and releases. 

Model: Our model consists of a square array of blocks, 
interconnected among nearest neighbors by coil springs 
with spring constant K and relaxed spring length I. To 
compare with experiments, free boundary conditions are 
used, with three nearest neighbors on edges and two 
at corners. Initially, the blocks are randomly displaced 
about their mean positions f* = (ia,ja) by (x, y) on a 
rough substrate, where i,j — 1, . . . ,L, and a is the lat- 



tice constant. Motivated by the drying process described 
above, internal stress may be introduced in various ways: 
a may be fixed but I is decreased in time to model con- 
traction; or both a and I may be fixed to impose an initial 
tensile strain (s = (a — I) /a > 0), while K is increased 
to model material stiffening Q . As a first study, we will 
adopt the simpler second approach here. 

Contrary to similar models of earthquakes |l^JTl[], in 
which the stress arises from relative motion of two sur- 
faces and released by block slips alone, the stress here 
is imposed uniformly by increasing K and released by 
block slips and spring breaks via a threshold dynamics: 
a block slips to a force-free position if the net force from 
its neighbors F > F s (the threshold for slipping), and 
a spring breaks (K = 0) if the tension b > F c = kF s 
(the threshold for cracking). With all forces initially be- 
low thresholds, K is increased slowly with F s and F c 
fixed. When either threshold is reached, stress is dissi- 
pated and redistributed accordingly, causing further slips 
and/or cracks until F < F s and b < F c everywhere. This 
constitutes one event. K is increased only between in- 
dividual events, corresponding to slow drying where the 
rate of driving is infinitesimal compared to that of relax- 
ation. 

The force exerted on a given block at f = f + (x, y) 
by a neighboring block at r' is given by (|r — r'\ — l)K. 
Since cracking is primarily due to tensile stresses, this 
nonlinear dependence on the coordinates leads to unnec- 
essary complications in updating the configurations. To 
simplify and compare with previous models, we expand 
the forces to first order in (x, y) to obtain the force com- 
ponents on a block at (cf. [pd| ) 

F x = (a - I + Xi+ij - Xij)Ki + - x l j)sK 2 + 

(-a + I + Xj-ij - Xi,j)K% + (xi t j-i ~ x l:j )sK4 

F y = (yi+i,j - yij)sKi + (a-l + yij+i - yi,j)K 2 + 

(Vi-ij - Vi,j) sK 3 + {-a + 1 + Vi,j-i - Vi,j) K i, (!) 

where the subscript of K denotes different springs. Like- 
wise, the tension of the spring between two blocks at (i, j) 
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and (i + has the components b x = (a — I + Xi+ij — 
x it j)K and b y = - Vi,j)sK. 

Without loss of generality, we hereafter choose our 
units such that a = 1 = F s . Of the remaining parameters 
|L2| {s, k, K}, K may be used to define the 'time' t = K. 
Note that a larger k means either a stronger material 
(larger F c ) or a weaker substrate coupling (smaller F s ). 

Results: We have simulated the model for a wide range 
of values of the dimensionless parameters: 0.01 < s < 
0.98, 1 < k < 30, and system sizes 20 < L < 300. 
Clearly, in the slippery limit k = oo, no spring is broken 
and our model becomes a variant of previous slip-stick 



models |H|,10| in their conservative limits, with almost 
the same slip-size statistics. However, notice that our 
model is driven multiplicatively by K, whereas previous 
models are driven by an additive force term. For finite 
k, the system contracts by block slips before it cracks 
into disjoint fragments. In agreement with intuition and 
experiments both the waiting time t c for the first 

spring to break and the mean fragment size increases with 
k (see Fig. [I] and below). 



(a) {0.5,1.6} 



(b) {0.5,2.3} 



with a long plateau and a sharp roll-off (see Fig. ||(a)), 
corresponding to a peak there in the probability den- 
sity D{c). The cutoff c max scales as L 2 . The distribu- 
tion of fragments is inhomogeneous deep in this region, 
as shown in Fig. ^(c). Along the phase boundary at 
k = k*(s), D > (c) ~ c~ r ' (see Fig. f§(b)) with an expo- 
nent rj S3 0.75 ± 0.05 to 0.63 ± 0.02 as s varies from 0.8 
to 0.5. 
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event. 

ft* S3 2.4; (b) power-law decay -D>(c) ~ c -0 ' 63 at ft = k 
(c) the phase diagram with a possible critical end point 



2. (a) Cumulative distribution of crack size per 
Note short (long)-range correlation below (above) 



FIG. 1. Typical patterns for L = 100. Legends stand for 
{s, ft}, (a) Below transition, shows localized cracks and small 
fragments; (b) near transition, almost percolating cracks and 
larger fragments; (c) inhomogeneity well above transition; (d) 
in the one-phase region, shows more isotropic cracks. 

Fig. [Ijalso reveals two kinds of crack patterns, depend- 
ing on an s-dependent ratio k* for large s: 

(a) Static cracking — For k < k* , cracks propagate slowly. 
Correlations in positions and orientations are weak be- 
tween successive cracks. The cumulative distribution 
D>(c) of the number of cracks per event, c, decays ex- 
ponentially, it is short-range (see Fig. ||(a)). The spatial 
distribution of fragments is homogeneous, as shown in 

Fig. g. 

(b) Dynamic cracking — For ft > k*, cracks percolate the 
sample and break it in a catastrophe. D > (c) is long-range 



By probing the k and L dependence of c max , we de- 
termine the phase diagram, Fig. ||(c). The transition 
between static and dynamic cracking becomes less strik- 
ing as s is decreased. Beyond the end point (s c , k c ) 
|13| , the two phases become indistinguishable. Travers- 
ing along fixed s within the one-phase region, we find 
in c max versus k a mild peak which saturates for large 
L, and D > (c) eventually decays exponentially at large c, 
indicating non-percolating cracks. 

Also of particular interest is the probability density 
P(f) of fragment area /, the center of focus of fragmen- 
tation theories [ 
percolating crac 
normal distribution ItLI 



14|,|15| and experiments [|16|. For non- 
ring, P(f) is best described by a log- 
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as shown in Fig. ^(a) for the cumulative distribution 
P> (/) = j^ 2 df P(f'). On the other hand, P> (/) for dy- 
namic cracking deviates from log- normal for k > n* . At 
larger n, a kink appears and becomes more pronounced 
(see Fig. ||(b)), implying different statistics for the frag- 
ments near the boundary and the center. Although the 
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tails of P> (/) suggest a power-law dependence, we can- 
not be certain due to the limited range of data. 
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FIG. 3. Cumulative distribution of fragment area, 
(a) shows good fits (lines) to log-normal distribution for 
non-percolating cracking; and (b) shows deviation from 
log-normal and emergence of a kink (spatial inhomogeneity) 
for dynamic cracking, where K = 2, 3, 3.5, 4, 4.5, 5, 5.5, and 
6 from left to right. 

A consistent physical picture emerges from study ing 
the reduced stress field e(r,t) = \b\/F c oc \/Tr cr 2 / F c , 
where cr is the local stress tensor, and < e < 1. We 
compute the structure factor S(k,t) = (\e(k, t)\ 2 ) / L 2 — 
L 2 d^ Q (e(r,t)} 2 , where e is the Fourier transform, and 
the over line means a spatial average. Averaging over 
the circular average S c i r (k,t), we define a characteristic 
length of the stress field R(t) = (k 2 )^^ 2 . Before crack- 
ing occurs, block slips start from the edges and propa- 
gate inward. Thus, e is reduced near the edges but en- 
hanced linearly in time in the bulk. These changes are 
reflected in the growth of R(t) via a power law (cf. j^]) 
R(t) <~ f . Within statistical uncertainty, <fi ~ 1 f° r aU 
s asymptotically at large L. The self-similar build-up of 
correlation is evident in the temporal scaling for t < t c : 
S cil (k,t) = ^$(fc^) (see Fig. |(a)). Hence, R(t) may be 
interpreted either as the penetration depth of slip events 
into the bulk or the correlation length of the stress field, 
whose growth towards L is a manifestation of the ap- 
proach to a critical state self-organized by the stick-slip 
mechanism EflO. 

The time t c when cracking sets in may be deduced by 
subjecting the enhanced stress |6| « Ks in the bulk to 
the threshold condition e « Ks/F c = ts/n > 1. Thus 
t c « k/s, which has been verified. Whether the initial 
crack triggers an instability depends on the strength of 
the stress S at a crack tip relative to the spatial fluc- 



tuation Se. Since Sb ~ KA, we have Se/e oc A/s. So, 
for small s, S is suppressed relative to Se, fluctuations 
dominate and the cracks are localized. This explains the 
absence of percolating cracks in the small-s regime. 

For large s, the transition can be understood by virtue 
of the effect of k on S. S is large only if there are stress 
transfers from crack sides to tips by slipping those blocks 
that link to broken springs. The ratio F/F s ~ Ks « k 
determines if they slip, where F ~ Ks in the bulk (see 
Eq. [l]) and K t c have been used. Since n has little 
influence on Se, we conclude that small k gives rise to 
small S and hence isolated cracks. The cracks border- 
ing a typical fragment are highly localized and weakly 
correlated, the usual argument leading to a log-normal 
distribution applies fl~^ , |l4|| . On the other hand, large k 
fulfills the slipping criterion easily, we have substantial 
stress transfers at a crack opening and a large S, which 
forces more broken springs. An instability is inevitable. 
The strong correlations of cracks (as shown in D > (c)) 
invalidates the assumption for log-normal distribution. 

Further insight can be gained from exploiting R(t). 
For non-percolating cracking, we find the mean fragment 
area (/) ~ R 2 , where R(t c (n)) = R c (k) (see Fig. f|(b)). 
This implies R{t) 2 prior to cracking sets a lower bound 
of (/) for the final cracked state. Using t c ss k/s, we find 
(/) oc k 2 ^. This simple relation highlights the crucial role 
of k in our model. In contrast, for dynamic cracking, this 
predictability is lost as e is extensively modified by the 
block slips at crack openings. 
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FIG. 4. (a) Scaling of structure factor for stress field prior 
to cracking, corresponding R(t) = 30-180. (b) Mean fragment 
area scales with correlation length of stress field for static 
cracking (upper curve), but not for dynamic cracking (high 
end of lower curve) . 

Finally, our model is robust in the following sense. If 
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we allow F c and F s to increase with K, the simplest al- 
ternatives are: (a) F c oc K, F s = const, so n oc K; (b) 
F s oc K, F c — const, so n oc 1/K; and (c) F s oc F c oc i'C, 
k = const. In the first case, the final state is entirely 
controlled by F c /K. It undergoes a first order transition 
from a crack- free to a cracked state as F c /K — > s+2A = £ 
from above. Further below the density of broken 
springs quickly saturates. The second case is trivial, hav- 
ing all the springs broken at long times. For the third 
case, either nothing happens if neither thresholds are ex- 
ceeded, or the updating sequence is ambiguous if both 
are. Thus, within our formalism, our choice of rules is 
the only one that evolves into a nontrivial cracked state 
through the interplay of sticking, slipping and cracking, 
without the need of fine tuning the parameters or initial 
conditions. 

Conclusion: We have introduced a robust spring-block 
model of surface fracture. Despite its simplicity, it cap- 
tures the correct trends observed in experiments and dis- 
plays rather complex behavior, which can be understood 
from the evolution of the stress field and characterized 
by dimensionless control parameters. Concerning the 
fragment distributions, similar transitions between log- 
normal and power-law distribution have been observed 
in some recent fragmentation experiments |n| . While it 
is possible that the underlying mechanism for the transi- 
tion may be similar to ours, one must be cautious because 
the physics with and without a substrate could be utterly 
different. 
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